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Abstract 



We develop a quantum kinetic theory of the chiral condensate and meson quasi- 
I particle excitations using the 0{N) linear sigma model which describe the chiral 

■ phase transition both in and out of equilibrium in a unified way. A mean field ap- 

proximation is formulated in the presence of mesonic quasi-particle excitations which 
^ \ are described by generalized Wigner functions. It is shown that in equilibrium our 

^ ■ kinetic equations reduce to the gap equations which determine the equilibrium con- 



densate amplitude and the effective masses of the quasi-particle excitations, while 
linearization of transport equations, near such equilibrium, determine the disper- 
sion relations of the collective mesonic excitations at finite temperatures. Although 
all mass parameters for the meson excitations become at finite temperature, ap- 
parently violating the Goldstone theorem, the missing Nambu-Goldstone modes are 
retrieved in the collective excitations of the system as three degenerate phonon-like 
modes in the symmetry-broken phase. We show that the temperature dependence of 
the pole masses of the collective pion excitations has non-analytic kink behavior at 
the threshold of the quasi-particle excitations in the presence of explicit symmetry 
breaking interaction. 
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1 Introduction 



In our previous paper [1], hereafter referred to as I, we have formulated a 
kinetic theory of self-interacting meson fields with an aim to describe the 
freezeout stage of the space-time evolution of matter in relativistic heavy- 
ion collisions. By using a single component real scalar field model we have 
obtained a set of coupled equations of motion for the meson condensate and 
the quasi-particle excitations. The meson condensate obeys the classical field 
equation with a modification due to the coupling to mesonic quasi-particle 
excitations, which is expressed in terms of the Wigner functions defined by 
quantum statistical expectation values of the bilinear forms of the quantized 
field operators. The Wigner functions contain the off-diagonal components 
which can be eliminated by Bogoliubov transformation in uniform systems, 
but survive in general non-uniform systems. We also solved the kinetic theory 
in and near equilibrium and studied a gap equation and the dispersion relation 
of the collective excitations. 

It is the purpose of the present paper to apply this formalism to the 0{N) 
symmetric linear sigma model which possesses a continuous symmetry. For 
= 4, this model may be considered as an effective field theory of QCD 
with light quark fiavors: the chiral condensate is described in terms of the 
sigma meson condensate, which is defined by the expectation value of the first 
component of the A'"-component scalar quantum fields. 

The new aspect which arises due to the presence of continuous symmetry is 
the manifestation of the Nambu-Goldstone modes[2, 3] which are identified as 
pions. The chiral phase transition has been studied using various effective chi- 
ral models of QCD, such as the linear sigma model[4] with flavor extension[5], 
the Nambu-Jona-Lasinio model for effective quark fields[6], and the non-linear 
realization of Weinberg as a low energy effecive theory of QCD [7]. Here we 
adopt the 0{N) sigma model which has been originally used in the context of 
electro- weak phase transition [8]. The model has been used by many authors 
to describe the chiral phase transition at finite temperature in the mean field 
(Hartrec) approximations [4, 9, 10, 11, 13, 15]. We shall extend such descrip- 
tions to non-equiibrium situations with the formalism developed in I. 

In this paper we first derive a coupled kinetic equations for the chiral con- 
densate and the mesonic quasi-particle excitation and show these equations 
reduce in equilibrium to the gap equations to determine the masses of the 
quasi-particle excitations in the Hartree approximation. The solution of these 
gap equations exhibits a first order phase transition. This is a well-known 
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generic problem of the mean field approximation which may be remedied by 
the inclusion of higher order fluctuations[12]. It has been known also, however, 
that such singular behavior of the chiral phase transition is smoothed out also 
by the introduction of the explicit symmetry breaking so that it generates the 
observed pion mass in the vacuum. In this work, we also study the effects of 
the explicit symmetry breaking in the present formalism. 

Another difficulty associated with the mean field approximation is that the 
pion quasi-particle excitation becomes massive apparently violating the Gold- 
stone theorem [3]. This is another generic problem of the Hartree approxi- 
mation. It has been known that the missing Nambu-Goldstone modes arise 
as acoustic collective modes of the excitations [20, 21, 22, 23]. We demon- 
strate this by computing the dispersion relations of the collective excitations of 
the system near equilibrium from the linearized kinetic equation with respect 
small deviations of the condensate amplitude and the quasi-particle distribu- 
tions from their equilibrium values These modes appears in the continuum 
of the quasi-particle excitations and thus may suffer a collisionless damping 
(Landau damping). We examine how these Nambu-Goldstoned modes are af- 
fected by the inclusion of the explicit chiral symmetry breaking. It will be 
shown that the Nambu-Goldstone mode becomes a massive excitation in the 
quasi-particle continuum which turn into un-damped pionic excitation in the 
high temperature symmetry recovered phase through a kink non-analytically. 

In the next section we extend the formalism developed in 1 to derive the 
coupled kinetic equations for N component self-interacting fields with 0{N) 
symmetry. We will show that N component meson condensate is described 
by N non-linear Klein-Gordon equations which contain non-linear interaction 
terms given in terms of the classical mean fields and N x N component ma- 
trix of fiuctuating fields. The quasi-particle degrees of freedom is expressed 
in terms of a 2 x 2 matrix of generalized Wigner functions which consists of the 
ordinary one-body density matrices in diagonal components and "anomalous" 
one-body density matrices in off-diagonal components, each component being 
N X N matrix of A^ component fields. For the long wavelength excitations, 
the equations of motion of the diagonal components of the generalized Wigner 
functions are shown to reduce to the Landau- Vlasov-type kinetic equations for 
the quasi-particle distribution function with a source/sink term expressed in 
terms of the off-diagonal components of the generalized Wigner functions. Ki- 
netic equations similar to ours have been derived for dilute cold atomic gases 
in [19]. 

In section 4, uniform static solutions of these coupled kinetic equations are 
shown to reproduce the well-known Hartree equilibrium solutions at finite tem- 
perature which exhibit the first order chiral phase transition. All off-diagonal 
components of the generalized Wigner function vanish in equilibrium and the 
A^ diagonal components give thermal distributions of the A" kinds of quasi- 
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particles, one being sigma-like and other N — 1 pion-like modes. 

In section 5, we apply our coupled kinetic equations to compute the disper- 
sion relations of the collective excitations near equilibrium. By linearizing the 
kinetic equations with respect to a small deviations from the static, uniform 
equilibrium solution we derive a set of self-consistency conditions to deter- 
mine the dispersion relation of the collective mode. It will be shown that the 
fluctuation of the diagonal components oi N x N matrix are all coupled and 
give dispersion relation of sigma-like collective mode while pion-like collective 
modes are contained in the off-diagonal fluctuations of the Wigner functions. 
The missing Nambu-Goldstonc modes are retrieved in these pion-like collec- 
tive modes as acoustic modes by making use of the self-consistency conditions 
to determine the condensate in equilibrium. We examine how these collective 
modes are modified when the chiral symmetry is not exact at the end of this 
section. 

A short summary of this work is given in the last section together with some 
open problems. The explicit expression of our kinetic equations in 0{N) model 
are given in the appendix. 



2 Kinetic equations for 0{N) linear sigma model 



In the previous paper [1], wc discussed coupled kinetic equations for one- 
component real scalar field, which has no continuous symmetry. In order to 
apply for more realistic physical situation, we need to extend the model to 
0{N) symmetric linear sigma model which possess continuous symmetry. The 
new aspects arise due to the presence of continuous symmetry, namely the 
manifestation of the Nambu-Goldstone modes, can be seen in the 0{N) sigma 
model. 



2.1 Heisenberg equations of motion for 0{N) sigma model 



The Hamiltonian of 0{N) sigma model is given by 



H 



dr 



(2.1) 



where scalar fields (pi and their canonical conjugate momentum field tTj are 
quantized by the usual equal-time commutation relations: 
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4>^{r, t), 4>j{r', t)] = [7r,(r, t), 7r,(r', t)] = 0. 



(2.2) 
(2.3) 



In the absence of the external field hi, H is symmetric with respect to rotaions 
of (f)i and Hi which form an 0{N) group. This Hamiltonian may be considered 
as a meson sector of the effective Hamiltonian of QCD if we set hi = e6i^i and 
e is chosen to reproduce the pion mass. 0i corresponds to the sigma meson 
field and (02, 03, 04) form isovector pion fields. 

The Heisenberg equation of motion of the quantum fields 4>i is given by 



dt 



— 1 



cl)i,H\=n,{r,t) (2.4) 
and the equation of motion of the canonical conjugate field 7r(r, t)i becomes 



N 



-i [tt,, H] = (V^ - m')Ur, t) - ^0,(r, t) 0?(r, t) - h, (2.5) 



dt 



After eliminating the momentum field vr from these equations, we obtained a 
Klein-Gordon equation for the quantum scalar field (pf. 



□0,(r, t) + m^4>^{r, t) = -^0i(r, t) 0|(r, t) - hi 



N 



(2.6) 



2.2 The mean field approximation and Gaussian Ansatz for fluctuations 



We define the mean field by the quantum statistical average of the quantum 
fields and then introduce as in I, the Gaussian Ansatz for the density operator 
p with respect to fiuctuations 



0c,i(r, t) = {4>i{r, t)) = tr(p0i(r, t)), 

7^c,i{^,t) = {7Ci{r,t)). 

0i(r, i) and 7ti{r,t) are given by 



(2.7) 
(2.8) 



0i(r,t) = 0i(r,t) - 0c,i(r,t) 
7^i{r,t) = 7ri(r,t) -7rc,i(r,t) 



(2.9) 
(2.10) 
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(0,(r,i)) = O 
(0-(r,i)> = O 



(2.11) 
(2.12) 



and 



N 



N 



+ Y: (2(0.(r, i)0,(r, t)f + {4>l{v, t)) {4>]{v, t 



i=l 



(2.13) 



This Ansatz for the density operator imphes a non-equiUbrium generahzation 
of the Hartree approximation in equihbrium. These shifted field operators obey 
the same equal-time commutation relations as the original field: 



0i(r,t),7rj(r',t) =i5ij5(r-r') 
0i(r,i),0,-(r',d =[fr,(r,i),fr,-(r',i)] = 



(2.14) 
(2.15) 



2.3 The Wigner functions 



To derive a kinetic equation for quantum fluctuation 0j, we first define the 
particle creation and annihilation operators as 



(r,t) = ^e^P- 



7fi(r,t) =i^e 



4p(^) - «i.-p(^) 



with 



The operators obey the usual equal-time commutation relations: 



(2.16) 
(2.17) 



(2.18) 



«i,p(^)'4p'W 



(2.19) 
(2.20) 



We introduce the Wigner functions of 0{N) model by 



6 



Fij{p, k, t) = (4p+k/2(^)«j,P-k/2(i)) 

Fij{P: k, t) = (ai,_p_k/2(^)4-p+k/2(*)) 

GijiP, k, t) = (4-p-k/2 W4p-k/2 W) 

Gij{p, k, t) = (flj -p-k/2(t)aj,p_k/2(^)) 



(2.21) 

(2.22) 

(2.23) 
(2.24) 



These four Wigner functions are not independent but are related each other 
by the commutation relations (2.20), e. g. 



Fij{p, k, t) = Fji{-p, k, t) + Sijdk^o 
Other relations are given in Appendix A. 
If we introduce a two-component notation of the operators 



A,p(t) = (a.p(t), aU(t)) and Aljt) = i ""''"^^^ 



>t it) ' 

the Wigner functions may be written in a matrix form as 



(2.25) 



(2.26) 



W,, (p, k, t) = (4p+k/2(^)^.-,P-k/2) 



Fij{p,k,t) Gij{p,k,t) 



yGij{p,k,t) Fij{p,k,t) 
The Fourier transforms of the Wingner functions are defined by 



(2.27) 



w,,-(p,r,t) = ^e-*--W,,(p,k,i) 



where (2.25) imphes 



^ fij{P,r,t) gij{p,r,t)^ 
^gij{p,r,t) fij{p,r,t) j 



(2.28) 



fij{p, r, t) = fji{-p, r, t) + 6i 



(2.29) 



As we have noted in 1, we consider in (2.18) as free parameters which may 
be chosen as different from the mass parameter m in the original Hamiltonian 
of 0{N) model. Generally, physical particle masses for interacting fields may 
differ from the mass parameters in the Hamiltonian due to the various kinds 
of dynamical effects. Here we focus on time-evolution of the nonequilibrium 
system where the physical particle mass cannot be defined globally to take a 
fixed value. Since a different choice of the mass parameter provides us with 



7 



a different definition for "particle excitations" of a nonequilibrium system, 
we sliould make proper use of the Wigner functions with an appropriately 
chosen mass parameter /Xj when we compare the results with observed particle 
distribution. 

Possible instability of the system is signified by the appearance of negative 
value of mass square if we determine it self-consistcntly. However, if we fix 
the mass parameters /^j, this instability, if exists, would manifest itself as an 
instability of the solution of the equations of motion for the condensate. 

Suppose we choose a different particle "mass" /i^ to define the creation and 
annihilation operators: 



where 



(2.30) 



,«pr_ 



7f,(r,i)=z^e*P 



<pW+<-pW 



«?p(^) - < -p(^) 



(2.31) 
(2.32) 



with 



4,p = Vp' + /^?- 

These new operators should also obey the commutation relations 



(2.33) 



<p(t),a2p,(t)]=4,5p,p' (2.34) 
'<,pit), <p'(^)] = [<p(^), 4At)] = 0, (2.35) 
and they are related to the original ones by the Bogoliubov transformation: 



(2.36) 



where the real parameter ctj^p is given by 



and Ti,T2, T3 are the Pauli matrix 
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,T-2 = 




,T-3 = 














,0-1, 



(2.38) 



The relations between the new Wigner functions defined by the new creation 
and annihilation operators {a^,ai) and the original Wigner functions are de- 
scribed as the following matrix equation: 

W^(p, k, t) = (4U.(t)4p„|(t)) (2.39) 
= /-.P+§ K p+t W^.,P-t(^))e"" ""^ (2-40) 
= e"'"''P+^ Wi,(p, k, 1)6""^'^-^ (2.41) 

For small deviations of the mass parameters /li — > /li + S/ii and /ij — > /ij + Siij, 
the change of the Wigner functions will be written by 



SWij{p, k, t)^TiWij{p, k, t)5ai,p+k/2 + Wij(p, k, t)ri5aj-p_k/2 
^ TiWjj(p,k, t)/i^(5//^ Wjj(p,k, t)nnj5nj 

2((p + k/2)2 + ;,f) ^ 2((p - k/2)2 + /.f ) ^ • ^ 



^.^ Equation of motion in the mean field approximation 



We obtain the equation of motion of the classical mean field (f)c,ii^, t) by taking 
the quantum statistical average of the Klein-Gordon equation of the scalar field 
of (f)i (2.6). With the Gaussian Ansatz for the density matrix, we obtain the 
equations of motion of the classical mean field 4>c,i {i = 1, 2, A'"): 



'n + m^)(j)cAr,t) = - 



A 
3! 



0. 



TV 

c,i X] 



+ (0?)) + 2 



N 



-hi (2.43) 



As we have done so in 1, we may call these classical equations non-linear 
Klein- Gordon equations of O(A^) sigma model corresponding to the non- linear 
Schodinger equation (or Gross-Pitaevskii equation) in the theory of the Bose- 
Einstein condensates. The non-linearity originates from the self-interaction of 
the classical fields (pc.i^ the interaction between 0c,i and (pcj, and also from the 
interaction with the fluctuations (0j(r, t)0j(r, t)) which implicitly depends on 
4>c,i- The fluctuations {(f)i{r,t)(f)j{T,t)) arise from "particle excitations" because 
of the following expression in terms of the Wigner functions as 
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(r,t)0,(r,t)) = ^ 
p>p' 

= E 



r,i{-p+p')-T 



:((aj,_p + 4p)(aj,p/ + al_p,)) 



p,p' V'^^iW^^3,P' 

gi(-p+p').rtj- [(1 + ^^)Wy (E±P:, p _ p', t) 



(2.44) 



Thus, the time-evolution of the classical mean field (pci is coupled with that 
of the Wigner functions. 



The equation of motion of the Wigner functions may be obtained from the 

commutato 
defined by 



commutators of bilinear forms of Oj^p and a] p with the mean field Hamiltonian 



Hrr,f= / dr 



i ij 



/4 



where 



" k=l 



2 [(pc,iYc,j 



and 



(2.45) 



(2.46) 



which implies that the mass parameters are given by 



(2.47) 



" i=i 



(2.48) 



The momentum representation of this mean field Hamiltonian can be written 

as 



^mf = E ^i.P^lp^i.P + O E ^n^.q 



p+q^ 



p,i 



p,q>2>j 



V^'^i,pV^'^j,P+q 



(2.49) 



where 
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An 



(2.50) 



In the Heisenberg picture, time-derivative of operators are calculated from the 
commutator of the Hamiltonian. With the Gaussian Ansatz for the density 
operator, the average of the commutator with the original Hamiltonian is 
equivalent to that with the mean field Hamiltonian H^f. 



(2.51) 



Therefore, we obtain the equation of motion of the Wigner function Wij(p, k, t). 
These equations can be written in a matrix equation as: 



d 

«^Wij(p,k,t) 



-^i,p+k/2T3Wij(p, k, t) + CJj-,p-k/2Wjj(p, k, t)T3 

T3(l + ri)Wi,(p + q/2,k + q,t) 



+E 



2^j,p+k/2 y 2c<Ji^p+k/2+q 

W,Kp-q/2,k + q,t)(l + ri)r3 



z,q ^2a;/,p_k/2 -y/2a;j,p_k/2-q 



AH 



(2.52) 



Eq. (2.52) and the nonlinear Klein-Gordon equations (2.43) form a closed sys- 
tem of coupled differential equations. We note that the effect of the exthernal 
field hi, only appears in modification of the nonliear Klein-Gordon equations. 



3 Kinetic equations for slowly varying system: the Vlasov equa- 
tions of 0{N) model 



In this section we derive kinetic equation of 0{N) model in the long wave- 
length limit. Here we assume inhomogeneity of the system is due to the long 
wavelength fluctuation k,q,<^p. We use the following approximations: 



p ■ k 

^i,p±k/2 ^ i^i,p ± (3.1) 

, p ■ k / 1 1 \ , , 

uJi,P+k/2 ± i^j,p-k/2 ^ {uJi,p ± c^i,p) + — T — (3.2) 



and 
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1 ^^_fi^£Ji^^it±M- 

We also use the Taylor series expansion of the Wigner functions: 



(3.3) 



W,,(p±, q/2, k + q, i) ~ W(p, k + q, i) ± -q • VpW,,-(p, k + q, t) (3.4) 
Then the equation of motion of the Wigner function Wij(p, r, t) becomes 



Wy (p, r, t) = i(a;p,iT3Wy (p, r, t) - u;pjWij-(p, r, t)Tz) 



1 



1 



+ ^ o • V(r3Wij(p, r, t) - Wij(p, r, t)r3) 



+r3(l + Ti)^ iC/ii(r,t)w2j(p,r,t) + 



4 la;/ 



+ 



hP ^i,P 



p-V{Uu{r,t)-Wij{p,r,t)) 



Hp 



{VUu{r, t))wij{p, r, t) + {VUa{r, t)) ■ VpWi,-(p, r, t) 



+ - iUij{r, t)wii{p, r, i) + 7 ( + ] P ' "^(Uijir, t)wii{p, r, t)) 



4 V^'p ^|,p. 



Hp 



{VUij{r, t))wuip, r, i) + (Vt/,,(r, t)) • VpWa(p, r, t) 



+ (3.5) 



where we have defined the generalized mean field potential Uij as 



Uij{r,t) 



An,,(r,t) 



In the long wavelength approximation, we also have 



(3.6) 



(0i(r,i)0,(r,i))~^ 



tr [(1 + ri)wy(p,r,t)] 



2^ 



(3.7) 



which appears in both the self energy term njj,(2.46) and the generalized 
mean field potential Uij (3.6) and in the non-linear Klein-Gordon equation for 
the condensate 4>c,i (2.43). 

These equations constitute the generalization of the Vlasov equation we de- 
rived for one component scalar field in I. The equations of each components 
of Wij are shown explicitly in Appendix. They have rather complex struc- 
ture. The diagonal components {i = j), of these equation can be rewritten 
in the form similar to the classical Vlasov equation. In fact, we extract the 
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terms which contain the diagonal components of the Wigner function. Writing 
fa = fi, Uii = Ui,we obtain 



/i(p,r,t) 



-d- ^^^1 P ■ Vr/,(p, r, t) + VrUiir, t) ■ Vp/,(p, r, t) 



+If + L + I-a 



'f^^g^^g (3-8) 
where If contains the terms proportional to the off-diaginal components of fij 



i \Uij(r, t)fji(p, r, t) - fij(p, r, t)Uji(r, t) 
1 



1 / 1 
+- 



4 V^'p ^Ip, 
1/1 1 ^ 



p • VrC/ij(r, t)fji{p, r,t)+p- Vr/ij(p, r, t)Uji{r, t) 



,\ 2 ' 2 I \Uij{r,t)p-Vrfji{p,r,t) +p-Vrfij(p,r,t)Uji{r,t) 

4 \^i,p ^j,pj \ , 



+ 2 ( VrC/i,(r, t) ■ Vp/,i(p, r, t) + Vp/^ (p, r, t) ■ V,U^,{y, t) 



while Ig and Ig are the terms linear in gij and gij, respectively: 



(3.9) 



iUij{r,t)gji{p,r,t) 



1/1 



4 V'^'p '^Ip. 



p ■ VrUij{r,t)gji{p,r,t) 



■7(^ + ^1 Uij{r,t)p-\/rgji{p,r,t) 

4 \'^i,p ^j,p/ 



--VrUij(r, t) ■ Vpgjiip, r, t) 



(3.10) 



-WijiP, r, t)t/jj(r, t) + i ( ^ J- ) 5ij(p, r, t)p • VrC/ii(r, t) 

+ 7 ( + ^ J P- Vr^ij(p,r,t) • Uji{r,t) 
4 V'^i.p ^j,p/ 



+ 2^p^^i(P' ^' ^) ■ '^r^^ii(r, t) 



Introducing the quasi-particle energy of the i-th species by 

ei{p,r,t) = uji^p + Ui{r,t) 



(3.11) 



(3.12) 
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the first two terms on the right hand side of (3.8) takes a form of the Landau 
kinetic equation: 



/i(p, r, t) = -Vp£j(p, r, t) ■ Vr/i(p, r, t) + Vr£i(p, r, t) ■ Vp/i(p, r, i) H 

(3.13) 

where the first term of the right is the quasi-particle drift term with drift 
velocity given by 

Vp£.(p, r,t)^ — (l- ^^M") p . Vr (3.14) 



The other equations for the time derivative of the off-diagonal components fij 
as well as gij and gij have no classical counter-parts. 



4 Statistical equilibrium and the gap equations 



In this section, we discuss equilibrium solution of coupled kinetic equation for 
0{N) model in spatially uniform system. 

Taking hi — eSu, we assume that only one component of the meson field 0i 
has non-vanishing expectation value in equilibrium: 



0c. = ^'^^^ = (4.1) 

Inserting this condition into the non- linear Klein-Gordon equations (2.43), we 
have 



^Veq = ~ (^0e\ + 2(0?)eq. + T^^M | 0eq " 6 . (4.2) 

We also assume 



(0l4)eq = (4.3) 

for J = 2, ■ ■ ■ , iV. In the chiral limit e = 0, (4.2) has non-vanishing solutions 
for (peqiy^ 0) oiily low temperatures which satisfy 
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171= — 



3! 



N 



eq. 



j/eq. 



(4.4) 



In order to ensure that the time derivative of the Wigner functions vanishes, 
we furthermore impose the condition that the off-diagonal components of the 
self-energy Alljj or the mean field potential Uij vanish in uniform equilibrium 



(4.5) 



and so do the all off-diagional components of the Winger functions: 



Finally we set the diagonal components of the mean field potential also vanish: 







(4.7) 



so that the mass parameter /Xj are chosen as 



(4.8) 



This last procedure eliminates the diagonal components of the gijj and giy 

Now the only non- vanishing components of the Wigner functions are diagonal 
components of jij 



fij — fij — feq.{^i,p)^ij 

for which we take Bose distribution function: 



(4.9) 



/eq.(<^i,p) — 



where /? = l/ksT is the inverse temperature. 



(4.10) 



The conditions (4.1), (4.5) and (4.8) imply that the mass parameters should 
be chosen as 



2 2 I ^ 

l^i =m + - 



N 



(l + 25,,i)0,\ + 2(0,2)eq. + E< 



(4.11) 



15 



where the thermal fluctuations of quantum flelds are given by the relation 
(2.44) as 



<^f>.. - E = / 1;=-^ = (4.12) 

where the dimensionless function I^^ (x) is given by [1] 



(2)/ 



_ /-^o k'^dk i 

^ !^ _ _ 1 2 In + f J - ^7) - + 0{x') (4.14) 



6 2 4 47r V8 4 7 327r 

where 7 = 0.57721 ■ ■ • is Euler's number. We call these equations the gap 
equations of the 0{N) model. 

Ahhough, by setting hi = eSu and choosing (4.1), we have committed breaking 
the original 0{N) symmetry in the direction of 0i field, we still have an 
0{N — 1) symmetry with respect to the rotions of (02, <p3, - ' ' , 4>n)- We still 
have degeneracy in the mass parameters 112 — fJ-s — • • • = A* at- This implies 
that the original gap equations consisting of N equations are reduced into the 
following two coupled equations: 



^il^m' + ^<Pl^ + ^{4>IU. + ^{N+1){^IU. (4.16) 
and the static Klein-Gordon equation (4.2) reduces to 



^''/'eq = -^(0eq + 3(^?)eq. + (A^-l)(02)eq.)</'eq-e . (4.17) 



The reduced gap equations (4.15) - (4.16) together with the reduced static 
Klein-Gordon equation (4.17) determine the values of /ii, ^2 and 0eq. at given 
temperature self-consistently. As we also noted in I, our gap equations are 
identical to the variational conditions to determine the mass parameters in the 
CJT composite operator effective potential in the one-loop approximation[10] 
if we ignore the renormalization effect due to the divergent vacuum polariza- 
tion which we have dropped out in evaluating the Wigner functions. 
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4-1 Exact chiral limit (e = 0) 



We first show numerical results of the temperature dependence of the solutions 
of the gap equations (4.15) - (4.16) in the absence of symmetry breaking 
external field (e = 0). The results are plotted in Fig. 1. 



A^lOO 




T/Tc [MeV] 

Fig. 1. Temperature dependence of the solutions of the gap equations for the scaled 
sigma mass fii/no and the scaled pion mass iJ,2 / fJ-o StS a function of of t = T/Tc. 

In the low temperature phase where 0cq. 7^ 0, the gap equations arc solved 
with the equation (4.4) obtained from the static Klein-Gordon equation or its 
reduced form 



(4.18) 



Since all the thermal fiuctuations (0^)eq. in the gap equations vanish at zero 
temperature, all mass parameters vanish except for /ii which is determined by 



^^l{T = 0) 



(4.19) 



where we have used the relation 



m 



(4.20) 



to determine the vacuum condensate 0o- In Fig- li th^ solution of the gap 
equations are scaled by the sigma meson mass /xq at T = for fii = /ii//xo 
and fi2 = /i2/Aio as a function oi t — T/Tc in the case oi N — A. 
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The two mass parameters /ii,/i2 and the condensate amphtude 0eq. vanish 
at "critical" temperature Tc given by 



72 



(4.21) 



{N + 2)X 



which is obtained by setting /xi = /X2 = 0eq. = in equations (4.15) and (4.16) 



In the high temperature phase, the condensate vanishes 0eq. = and aU the 
mass paraeter become degenerate with 



which may be obtained from either (4.15) or (4.16) by setting 0eq. = 0, /^i = 



The phase transition in this model becomes of first order showing a typical 
hysteresis behavior. This behavior has been known as a generic symptom of 
the mean field approximation with the 0{N) sigma models with 0*^ interac- 
tions [4, 9, 11, 13, 14]. It has been shown, however, that the transition becomes 
second order if one includes the effect of the fiuctuations in the two-loop ap- 
proximation [12] based on the optimal perturbation method developed in [11]. 
The second order phase transition has been obtained also by the Nambu- 
Jona-Lasinio model[17, 18]. Although, in the chiral limit, the pion mass /i2 is 
expected to vanish due to the Goldstone theorem, 112 has non-vanishing val- 
ues in the low temperature phase. This apparent violation of the Goldstone 
theorem always occurs in the mean field approximation. With a resummation 
method known as an optimized perturbation theory at finite temperature de- 
veloped in 0{N) (j)^ theory, the Goldstone theorem is satisfied for arbitrary 
A^[ll]. In our framework, the missing Goldstone theorem can be retrieved as 
the collective excitation mode [20, 21, 22, 23] as we shall show in the next 
section. 



4.-2 With symmetry breaking (e 7^ 0) 

In the presence of the explicit symmetry breaking (e 7^ 0) the mass parameters 
IXi for i 7^ 1 also takes non-vanishing value even at zero temperature. This is 
seen by setting {(p'^)eq. = for all components of meson fields in (4.4) which 




(4.22) 




yields 
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+ \<i>l = ^ (4.23) 

where 0o is the amphtude of the vacuum condensate. Inserting this relation 
in (4.16) gives 



lA = j- (4.24) 

00 

at zero temperature. In the case of = 4 we may interpret this mass as the 
degenerate mass of three pions (7r+, ttq, 7r_) which form an isovector field. 
We thus set 



e = m^0o (4.25) 

as a condition to choose the symmetry breaking parameter e. One can also 
show that the equation of motion (2.6) imply that there is a partially conserved 
"axial vector" current 



dSid^^i - kd^'k) = 4i (4.26) 

Sandwitching this operator relation by the vacuum state and the single charged 
pion state we obtain the relation which connect the vacuum expectation value 
of (T field 01, namely 0o, to the charged pion decay constant J^, 



00 = h (4.27) 
The (7-meson mass in the vacuum is determined by (4.15) at zero temperature 



ml^m' + \4>l = + \fl (4.28) 



We use three conditions (4.25), (4.27) and (4.28) to determine three parame- 
ters rn? , A, and e in our model Hamiltonian and solve the gap equations (4.15), 
(4.16) and (4.17) to determine the temperature dependence of /Xi, ^2 and 0eq.. 
Numerical values of the parameters are determined by these conditions to 
give m„ = = 500MeV, = ^2 = MOMcV, = (pQ = 94MeV at zero 
temperature by tuning A = 78.2. We show numerical result computed with 
these values of parameters is shown in Fig. 2. It is seen that the transition is 
smoothed out with no trace of the first order transition. 
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Fig. 2. Temperature dependence of the solutions of the gap equations with exphcit 
symmetry breaking. 

5 Dispersion relations of the collective excitations 

In this section we compute dispersion relation of the collective excitation in 
the system near equilibrium in 0{N) model. We have shown in I that the 
collective mode appears in the low temperature phase as a meson mode with 
small effective mass as a coupling of the mesonic excitation with the quasi- 
particle excitations. In the 0{N) sigma model there are different kind of 
quasi-particle excitations in the system in equilibrium. We will show that the 
fluctuation of the condensate amplitude 0c,i couple with the fluctuation of 
the diagonal components of the Wigncr function fa, ga and gu and generate 
the collective mode of the sigma meson type, while the fluctuations of the 
condensate perpendicular to the condensate 0c,i with i > 2 couple with the 
off-diagonal components of the Wigner functions and generate pionic collective 
modes, which constitute the missing Nambu-Goldstone modes in the chiral 
limit. 

We first assume that each of the meson field and the distribution functions 
consists of uniform equilibrium term and a small fluctuations around it: 



where we have used the equilibrium conditions: f^^' = /cq.(^i,p)^ii) 9if = 0- 
We further assume that these deviations are caused by a perturbation in the 
symmetry- breaking external fleld Shi{T,t). 





fij{p, r, t) = /eq.Kp)'5ii + Sfijip, r, t) = fji{-p, r, t) 
9ijiP, r, t) = 5gij{p, r, t) = g*i{p, r, t) 



(5.3) 



20 



Then, by linearizing the non-hnear Klein-Gordon equation (2.43) with respect 
to S(f)i^c{^,t), we obtain 



;n + /i?)<50e,i(r,t) 



A 
6 

A 
3' 



N 



30eq<5(0?(r,t))+0eqE'^(^|(r,^)) 



J=2 



+ Shi{r,t) 
(5.4) 



)eq5(</.i(r, t)<Pi{r, t)) + Shi (for i > 2) (5.5) 



where the deviations of the fluctuations of the quantum field (f>i{r,t) from its 
equilibrium value are given in terms of the sum of the four Wigner functions: 



5(0,(r,O0,(r,i)) = E 



1 



1 



+5gij{p,T,t) + 5gij{p,r,t)] 



+Sgij{p, r, t) + Sg*i{p, r, t)] (5.6) 

where in deriving the last line we have used the following symmetry of the 
Wigner functions: 



Sfij{p,T,t)=Sfji{-p,T,t) 

S9ij{P,r,t)^Sg*i{p,r,t) 



(5.7) 
(5.8) 



The space-time dependence of the Wigner functions are determined by the 
linearized equations of motion which read for Sfij and Sgij: 



11 1 

SfijiP, r, t)^i{ui,p - uj^pjSfij - -( + )p • (VSfij) 

(/eq.(<^i,p) - /eq.(<^j,p)) 



«.p i.p 



+ -iV(5C/ij • Vp(/eq.(a;i,p) + /eq.(t^j,p)) 



(5.9) 



11 1 

(^^ij(P, r, t) = -i(a;i,p + ^^^p)Sgij + - ( )p • (V^^^j) 



"^i.p "^j.p 

+ -V5C/i^- • Vp(/eq.(c^i,p) - /eq.(^^j,p)) 



(/eq.(Wi,p) + /eq.('^j,p)) 

(5.10) 
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and the equations of motion for 5fij and Sgij can be obtained from these 
by making use of the symmetry properties, Sfij{p,r,t) — Sfji{—p,r,t) and 
Sgij{p,r,t) = dg*i{p,r,t). 

The deviation of the mean-field potential SUij is caused by the change in the 
change in the condensate and the fiuctuation: 



5t/.,(p,r,t) = i^^M (5.11) 

with 



AT 
k=l 

+ ^ (0eq.(^0c,i(r, t)5ij + 0eq.50cj(r, t)5i^i) , (5.12) 

We solve the linearized equations of motions for a monochromatic perturbation 

5hi{r, t) = e-^'^+^+^-^^/ii + c.c, (5.13) 
which is adiabatically switched on (a;+ — uj + ie), with the following Ansatz: 

(50e,*(r, t) = e-*^+*+*'^-"(50i + c.c. (5.14) 
5/,,(p,r,t)=e— +*+^'^-5/,,(p) +C.C. (5.15) 
5gij{p, r, t) = e-^'^+*+^''--%,(p) + c.c. (5.16) 

where c. c. implies the complex conjugate. Resulting linearized kinetic equa- 
tions are shown in appendix D. From these equations, we find two kinds of 
dispersion relations of collective excitations near equilibrium. One of the re- 
lations is essentially equivalent to the dispersion relation of one-component 
scalar model obtained in I: it gives the dispersion relation of the sigma meson 
like collective mode. The other relation contains pion-like modes which can 
be interpreted as the Nambu-Goldstone modes of 0{N) linear sigma model in 
the chiral limit. 



5.1 Collective sigma modes in the chiral limit 



As seen in the linearized non-linear Kein-Gordon equation (5.4) the fluctuation 
of the sigma meson condensate (50c,i(r, t) caused by the external perturbation 
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5hi{r,t) couples with the change in the diagonal components of the fluctua- 
tions 5{4>j{r,t)). Along with the Ansatz (5.15) and (5.16) we put 



(5(02(r,t)) =^,e-^'^+*+^^-"- + c.c. 



(5.17) 



so that 



Prom (5.9) and (5.10) we obtain 



(5.18) 



/3v^,p • k (1 + /cq.(^»,p)) /eq.(^i,p) ,rT 

w+ - Vi,p ■ k 



where Vj^p = Vpa;j^p = p/<^i,p and we have used 



(5.19) 
(5.20) 



Vp/eqX'^i.p) = -Vi,p/3 (1 + /eq.(^i,p)) feq.{'^i,p) 



(5.21) 



while the change of the mean field potential is given by 



12u;,-, 



N 



20eq(2(5i,, + 1)501 + 2jr. + ^^_^. 



and the deviation of the sigma meson condensate is given from (5.4) 



(5.22) 



501 



6 



i=2 



(5.23) 



These equations form a closed set of equations to determine the linear response 
of the system to the external perturbation Shi. In finding the solutions it 
greatly simplifies to observe that the solutions have degeneracy due to the 
symmetry of the equilibrium solutions, 1^2 = • • • = fJ^N which implies J^2 ~ 
• • • = JT^v and El2-^i = {N- 1)^2- We find from (5.22) and (5.23) 



23 



5U^ 



A 



11 



12(x;i,p 



5U' 



A 



22 



12a;. 



2,P 



1 + 



(3^1 + (iV - 1)^2) 



hi 



+ 



A;2 



2 

■3' 



hi 



(5.24) 



(5.25) 



Inserting (5.19), (5.20) for i = 1,2 in the right hand side of (5.18) and using 
(5.22) and (5.23) we find self-consistency conditions for J^i and ^2'- 



^1--^ [/Ci(a;+, k) (3:^1 + (TV- 1)^2) + Ji(c^+,k)]$(u;+,k,//i) (5.26) 
J^2 = [ICi{uj+, k) (3j^i + {N- 1)^2) + 2{N + 2)^2 + Ji(w+, k)] 



x$(u;+,k, //2), 



(5.27) 



where 



/Ci(u;,k) = l + 
Ji(a;,k) 



(^2 _ Jj.2 _ ^2 
0eq.^l 



(^2 _ Ij.2 _ ^2 



(5.28) 
(5.29) 



may be considered as renormahzations of the quasi-particle interaction and 
the couphng to the external field, respectively, due to the couphng to the 
fluctuation of the meson condensate and ^{co^, k, /i) is the function introduced 
in I: 

$(cj+, k, fi) = $1(0;, k, n) + i$2('^, k, Ai) (5.30) 
with the real and the imaginary parts 



$i(a;,k,/i) = -P J 



d^p 1 



4c^a 



!(^,k, /^)=y" 



1 



(Vp • k)^ 



[27r (5 (a; - 2a;p) - 5{u + 2u^)) /eq.(^p) 



(5.31) 



+7r/?Vp • k {5{UJ - Vp • k) - 5{UI + Vp • k)) (1 + /eq.(a;p)) /eq.(c^p)] 

(5.32) 
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respectively, where we have indicated the quasiparticle mass // dependence 
exphcitly. 

Eqs. (5.26) - (5.27) forms two coupled linear equations of and T2 with an 
inhomogenious terms proportional to J^i which may be written in the matrix 
form: 



with 



M 



^1 
To 



Q2 



Ji 



(5.33) 



/ l + f/Ci$(l) 



12 



/Ci(Ar- 1)$(1) 



^/Ci$(2) 1 + A /Ci(A^-l) + 2(Ar + 2)U(2) 



12 



(5.34) 



and 



Qi 
Q2 



(5.35) 
(5.36) 



where we have used abbreviate notations $(1) 
$(a;,k, 112). 

Solving (5.33) for J^i and ^2, we find 



$(a;,k, /^i) and $(2) 



T2 



^(l + A(;V + 2)^Q,J, 



1 



detM 



Q2J1 



(5.37) 
(5.38) 



where 



detM= (^l + ^/Ci$(l)j (^1 + A(Ar + 2)$(2)j + ^^k:^{N - 

(5.39) 

These relations determine the linear response of the fluctuations 5(0^) to the 
external perturbation h\. Prom (5.23) 



501 = 



A 0, 



6 a;+ — A;^ — //f 



(3^1 + (TV - 1)^2) + .2 .2 ,.2 '^^i 



a;+ — A;^ — //^ 



(5.40) 
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which upon insertion of (5.37) and (5.38) gives 



501 = 



^5h 



(5.41) 



LO^ — — 111 



1 



with 



A0eq. [m (l + ^(iV + 2)$(2)) + ^$(2) 



(5.42) 



24{ujI - - 111) detM 



The quantity Si represents the mcdlTim modification of the propagation of the 
sigma meson field due to the absorption and the recmission of the fiuctuation 
by the pair of quasi-particle excitations. We note that Si is regular on the 
quasi-particle mass-shell cu'^ — k'^ + ix\ while new singularity appears when the 
condition 



is satisfied, corresponding to the collective excitations of the system. 

We have shown in I that the undamped collective mode appears in the time-like 
region in the low temperature phase where there exists non-vanishing sigma 
meson condensate (0eq. 7^ 0). To eliminate the singularity of ICi at meson poke 
we multiply (5.43) by uj'^ — k"^ — iif. One then finds at A; = 



for the condition of the shifted meson pole. 

The entire space-like region uj < k is covered by the quasiparticle continuum 
associated with the scattering of the thermally populated quasiparticle exci- 
tations which gives the non- vanishing imaginary parts when the kinematical 
condition cu — e^^p+k — e^^p ~ k • p is fulfilled. Since Vi^p — p/oJi^p ~ 1 — /Xj/p, 
this condition is met near the light cone uo = k only for quasiparticles with 
large values of p or cjj^p whose number is suppressed by the Boltzmann factor 
/eq.('^i,p) Indeed, we have found in I that the integral over p in (5.32) can be 
carried out analytically, yielding 




(5.43) 




(5.44) 
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Fig. 3. The region where the imaginary part of ^{u + ie, k) has non-vanishing value 
1 CO 1 



$2(w,k,/i) 



o_-2 U t^l^ 



— t(^{^ - 



(5.45) 



The first term gives continuum in the whole space-like region which diminishes 
exponentially as one approaches to the edge u = k due to the Bose-distribution 



factor /( 



eq. \^i,p) 



for ojr 



{uj/k)p = /i/y'l - {u/kf 



The second term in (5.45) may be interpreted as due to the thermally induced 
pair creation or annihilation, namely one of the pair of quasiparticles created 
or annihilated has been present as thermal excitations. It gives continuum in 
the time-like region oj = cUp+k -|- cUp ~ 2ujp > 2/i. This continuum thus covers 
the region for the uj greater than the twice the lightest quasiparticle mass, 
namely /i2 Since the sigma meson mass fii may be greater than two times the 
pion-like quasiparticle mass fi2, the collective mode of the sigma meson like 
mode has additional thermal decay width in addition to the pair decay width 
in the vacuum. 

For a long wavelength, low energy excitation mode (tu. A; <^ /ij), /Ci(u;, k) may 
be approximated as 



/Ci(cu,k)^/Ci(0,0) = l-— 1^ (5.46) 



In the exact chiral limit (e = 0), 
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- { J*^"- ^ °! (5.47) 



^0 (0eq. = 0) 

SO that we have 



^i(0,0) = <; (5.48) 



which gives 



det M = 1 - ^$(1) + ^$(2) - ^(^^ + 2)$(1)$(2) (5.49) 
for low temperature symmetry broken phase (0eq. 7^ 0) and 



det M = + A(iv + 2)$(1) j - ^$(1) j (5.50) 

for high temperature symmetric phase {(peq. = 0) where $(1) = ^(2). 

We show in Fig. 4 the numerical result of the shifted pole mass fi[ of the 
collective sigma-like excitation in the long wavelength limit {k = 0). At zero 
temperature ii[ coincides with /ii. As the temperature increases the coupling 
to the thermally excited quasiparticles makes the pole mass iJ,[ smaller than 
Hi . The mass shift exhibits non-analytic kink behaviors as it crosses the quasi- 
particlc threshold at c<j = 2/i2 and at certain temperature Ti slightly above Tc 
it becomes zero indicating the spinodal instability as discussed in 1. This insta- 
bility coincides with the condition for the appearance of the aforementioned 
tachyonic mode. 



5.2 Collective pionic modes in the chiral limit: the Nambu-Goldstone modes 



We now compute the collective excitations carrying the pionic quantum num- 
ber which couple to the oscillations of the meson condensate in the direction 
perpendicular to the equilibrium condensate, namely 50^. First we show that 
the missing Nambu-Goldstone modes appear as collective excitations with 
phonon-like dispersion relation in the exact chiral limit. 

The equations of motion of 6(f)i , (5.5), indicate that pionic modes couple to the 
fluctuation of the type 6{4>i4>i). There are N — 1 such mode for i = 2, . . . , A^, 
but these modes are all degenerate due to the unbroken 0{N — 1) symmetry, 
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Fig. 4. The spectrum of sigma-like excitations with k = 0. The effective mass fi'i of 
the collective sigma mode in the chiral hmit (e = 0) is plotted by thin solid curve. 
The continuum of the sigma-like quasi-particle excitations are shown by the gray 
areas for low temperature phase (T < Ti) and for high temperature phase (T < Ti) 
where Ti is highest temperature where the gap equation has non-zero solution for 



corresponding to the isospin symmetry in case of = 4. So it suffices to 
compute the i = 2 component to obtain the pionic excitation spectrum. 



Applying a monochromatic perturbation with non- vanishing 502, we set 



5(0i02)(r,t)=^i2e-*'^+*+*''-'- + c.c. 



(5.51) 



with 



12 



E 



2^t^l,pt^2,p 



^/i2(p) + 5/i2(p) + Sguip) + Sguip] 



(5.52) 



From the linearized equations of motion of the corresponding Fourier compo- 
nents of the Wigner functions given in Appendix D, we obtain 



<^/l2(p) 


=-^y(^+, 


5/l2(p) 




^^12(P) 




^^12(P) 


= -^12 (^+> 



,Kp)M2\p)SU^2 



(5.53) 
(5.54) 
(5.55) 
(5.56) 



where 
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(5.57) 



- ('^i.p =F cc;2,p) - k ■ (vi^p ± V2,p)/2 

A/-W(k, P) = (l + ^ - (/eq.K,p) T UM) 

+ (k/2)-Vp(/eq.Kp)±/eq.(^2,p)) 

(5.58) 

are the propagators of a pair of quasiparticles of two different types and the 
effective number of such states, respectively. The off-diagonal components of 
the disturvance caused in the mean-field potential may be obtained from (3.6) 



5C/i2(p) 



511 



12 



2^a;i,pa;2,p 



(5.59) 



with 



(^If 12 = - (0eq.^02 + ^12) 



(5.60) 



Using the solution of the linearized equation of motion for a classical conden- 
sate (5.5) 



X 



'■>(i(\.-r 12 



(5.61) 



we find 



(^C/12 = \ , ^ (/C2(u;+,k)^i2 - J2(c^+,k)) 



6 v^tc;i,ptc;2,p 



where 



X:2(a;,k) = l + - 



A </'e\. 



^2(0;, k) = 



3 u;2 _ k2 _ ^2 
1 



(5.62) 



(5.63) 
(5.64) 



Inserting (5.62) into (5.65) we find the self-consistency condition for the pionic 
excitations: 



:^i2 = Y^*(^+, k) [/C2(a;+, k)JF,2 - J^{u:^, k)] 



(5.65) 
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where 



1 



-A{uj,k, p) 



(5.66) 



with 



A{u, k, p) = {g{2\-u^, k, p) + gi^\u;, k, p)) Afi^\k, p) 

+ {g{t\-io, k, p) + g[t\io, k, p)) Mp{k, p) (5.67) 

is a modification of the dimensionless function $(u;,k) which has appeared 
in the computation of the sigma meson collective modes. Indeed, if we set 
1^1 = 1^2 = f^, the function \E'(co', k) coincides with ^{u, k). 
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Fig. 5. The region where the imaginary part of ^{uj + ie, k) has non-vanishing value 
is shown by gray area. 

Solving (5.65) for J-'i2 we find the expression for the linear response of the 
fluctuations to the exterrnal perturbation J^2'- 



^^(^+,k) 



12 



^^(^+,k)/C2(^+,k; 



-J2 



(5.68) 



The condition that Eq. (5.65) contains non-trivial solution gives the dispersion 
relation 



^/C2(^,k)vi/(^,k) = l 



(5.69) 
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Now we show that this equation has the solution of the form 



a; = v-k (5.70) 

which may be identified as the missing Nambu-Goldstone modes. 

For this purpose, we first show that the dispersion relation (5.69) is satisfied 
for (cj, |k|) = (0,0). We recall that In the chiral limit (e = 0) the static 
non-linear Klein-Gordon equation (4.2) and the definition of the quasiparticle 
masses (2.48) read 



m 



,.+3(0?)eq. + (iV-l)(02)eq. 



+ ^ (3<. + 3(0?)eq. + {N- l)(0^)eq.) 
^^' + ^«. + (0?)eq. + (iV+l)(0^)eq.) 



(5.71) 

(5.72) 
(5.73) 



respectively. Inserting (5.71) into (5.73), we find 



A 



((02)cq. - (0?)cq. 



(5.74) 



From (5.72) and (5.72), we also find 
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Fig. 6. The spectrum of pion-like excitations with A; = 0. The Nambu-Goldstone 
mode in the low temperature broken-symmetry phase at T < Ti is shown by the 
thick solid line as well as the massive pion mode degenerate with sigma-meson 
mode in the high temperature phase. The continuum of the sigma-like quasi-particle 
excitations are shown by the gray areas for low temperature phase (T < T\) and 
for high temperature phase (T < Ti). 
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^1-^1 = 1^ (0eq. + (<^?)eq. " (02)eq.) (5-75) 



Then, 



^2(0,0) = 1 - -—F = ,7ov 77:;^ (5-76) 



3 ^xi {4>i),^, - 



eq. 



where we have used (5.75). On the other hand, 



^(0, 0, p) = (6?i2 ^(0, 0, p) + g[:,\Q, o, p)) Ar^t^o, p) 

+ (0, 0, p) + g'^^ (0, 0, p)) Ulp (0, p) 

^ 2 /eq.(^l,p) -/eq.(^2,p) _ /eg. (^l,p) + /eg. (^2,p) ^ 
V <^l,p - ^2,p t^l,p + t^2,p 

^ 4(a;2,p/eg.(^l,p) - C^l,p/eg.(^2,p)) 

2 2 



(5.77) 



where we have used ^i?(0,0,p) = (cui^p =F cua.p) \ A/'/2^^(0,p) = /cq.(t^i,p) =F 
/eg.('^2,p) and p — a;|p = A*i — /^i- Inserting (5.77) into (5.66) we obtain 



^12 (0f)eq. - (0i)eq. 

A ■ + (0f)eq. - (0i)eq. 

where we have used (5.75) again. Comparing (5.76) and (5.78) we finally find 



(5.78) 



A*(0,0) = /C2(0,0)-^ (5.79) 

which imphes that the equation (5.69) has the solution (a;, k) = (0, 0). 

We note that the acoustic Nambu-Goldstone modes appear in the continuum 
of of the quasi-particle excitations where the function \I'(cJ+,k) gives non- 
vanishing imaginary part (see Fig. 5) thus always suffers the Landau-damping 
at finite temperature. 
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5.3 Collective modes with explicit symmetry breaking (e ^ 0) 

As we have seen in section 4, in the presence of the exphcit symmetry breaking 
e 7^ 0, the mass parameter /Xj becomes non-vanishing at zero temperature and 
we have interpreted it as the physical pion mass in case of = 4. At finite 
temperatures, the pionic collective modes, which form the masslcss Nambu- 
Goldstone modes in the chiral limit, thus acquire non- vanishing masses. These 
modes would also persist in the high temperature region since the transition 
becomes smoothed out as we have seen in the previous section. However, there 
may be some qualitative change in the character of the pionic collective modes 
in the low temperature region and in the high temperature region since the 
NG bosons appear in the chiral limit only in the low temperature phase. We 
now examine these problems. 

5.3.1 Collective sigma modes with symmetry breaking 

As shown in section 5.1, we have computed collective sigma modes with ex- 
plicit summetry breaking (e ^ 0). We show in Fig. 7 the numerical result 
of the shifted pole mass ^'^ of the collective sigma-like modes in the long 
wavelength limit. At zero temperature /x'^ is equal to the quasi-particle mass 
jii. As the temperature increased, shghtly exceeds //i at T < 120MeV, 
while it becomes smaller than /xi at T > 120MeV. The sigma mass shift has 
non-analysic kink similar to that in the case of e = when it crosses the 
quasiparticlc threshold at lo = 2/i2. However, it exhibits no instability since 
the first order transition is smoothed out. The shifted sigma pole mass also 
increases above T — 230MeV where the chiral symmetry is restored. 

5.3.2 Collective pionic modes with symmetry breaking 

We have computed the mass shift of the collective pionic excitations with ex- 
plicit symmetry breaking. By using Eq. (5.69), we show in Fig. 8 the numerical 
results of the pionic pole mass ji^ in the long wavelength limit. At zero temper- 
ature is equal to pionic quasi-particle mass ^2- As temperature increses 
becomes smaller than 112- At 150MeV< T < 170MeV, it exhibits non-trivial 
behavior. It has non-analytic kink where it crosses the quasiparticle threshold 

Recently, the shifted pole masses of the sigma/pionic collective modes has 
been investegated also by Tsue and Matsuda[24] based on the earher work[23] 
in the functional Schrodinger picture[25, 26, 27]. These authors discuss the 
cut-off dependence of the mass shift which arises due to the inclusion of the 
vacuum polarization, the effect which we have neglected in this paper. The 
kink structure of the mass shift is smoothed out in their results, however. This 
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Fig. 7. The spectrum of sigma-like excitations with explicit symmetry breaking. 
Effective mass ^'^ of cohective sigma meson mode with exphcit symmetry breaking 
is shown by the sohd Une. The gray area indicates the continuum of quasi-particle 
excitations with the same quantum number as a sigma meson. 
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Fig. 8. The spectrum of pion-Uke excitations with symmetry breaking. Effective 
mass ^2 of collective pion mode with explicit symmetry breaking is shown by the 
solid line. The gray area indicates the continuum of quasi-particle excitations with 
the same quantum number as a pion. 

may be due to the static approximation [uj = 0) they have taken in evaluating 
the pion self-energy which misses the threshold behavior of the quasi-particle 
excitations. 



To clarify this point, we examine the dispersion relations of the long wave- 
length pionic excitation modes for e 7^ 0, plotting in Fig. 9 
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Fig. 9. Graphical determination of the effective pion mass for the case of expUcit 
symmetry breaking. The function defined by Eq. (5.80) is plotted against uj for = 
at three different temperatures as indicated in each frame. The crossing points with 
the horizontal axis give the effectve masses of collective pionic excitations. 
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Diriuj.T) = lim Re 



Ax:2(a;+,k)*(a;+,k)-l 



(5.80) 



as the function of the excitation energy cu scaled by the temperature. We 
also plot in Fig. 10 the real part of the function ^{u;,k) which are used to 
evaluate as well as the imaginary part which determines the region of 
the quasi-particle continuum. The value of u at the crossing points of this 
curve with the horizontal axis give the energies of the collective excitations. 
At T — 128MeV the graph of D^^ has two kinks corresponding to the quasi- 
particle threshold //i ± /i2- As temperature increases the kinks shift to the left 
and the kink at fii — 112 exceeds the horizontal line. At T = 160MeV, therefore, 
we have three crossing points which corresponds to the non-trivial behavior of 
the pionic shifted mass as shown in Fig. 7 at 150MeV< T < 170MeV. At low 
temperatures, the tip of the kink is located below the horizontal axis. As the 
temperature exceeds, the tip of the kink moves upward and hit the horizontal 
axis giving two extra solutions of Dt^{uj.T) = 0. This is the reason why a pair 
of solutions appear with kink in the plot of the pion masses at the threshold 
of quasi-particle excitations aX. uj = fix — 112- 

We note that in the chiral limit the branch which originate from /i2 at T = 
becomes the massless Nambu-Goldstone modes. Our results show that with 
the explicit chiral symmetry breaking, the Nambu-Goldstone mode becomes 
a massive Landau-damping mode at low temperature which terminates at the 
kink position and there appears a new mode which becomes the pionic mode 
degenerate with the sigma meson mode in the high temperature symmetric 
phase. 



6 Concluding remarks 



In this paper we have developed a qiiantum kinetic theory for the chiral con- 
densate and quantum meson excitations based on the formalism we have de- 
veloped in our previous paper [1] applying to the 0{N) sigma model. We have 
shown that the time-evolution of system is described by a coupled form of 
kinetic equations consisting of the N classical non-linear Klein- Gordon field 
equations for classical meson fields (condensates) and 2N x 2N kinetic equa- 
tions for the generalized Wigner functions of quantized meson fields. 

We then applied our kinetic theory to describe uniform equilibrium state and 
to calculate the dispersion relations of collective modes near equilibrium. We 
constructed time- independent solutions of coupled equation for 0{N) model 
assuming only one component of the meson fields has non-vanishing expecta- 
tion value in equilibrium. We recover the results in the Hartree approximation 
at finite temperature. The order of the phase transition, calculated from the 
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Fig. 10. Real part (upper panel) and the imaginary part (lower panel) of the function 

•^{0;, k) at k = 



gap equations, becomes first order in the chiral limit, while it is smoothed out 
when we introduce symmetry breaking in order to generate actual pion mass. 
It is also known that the Goldstone theorem (in the chiral limit) is apparently 
violated because all mass parameters in quantized field fluctuations become 
non-vanishing even in the low temperature symmetry broken phase. 



It was shown that the dispersion relation of 0{N) model decouples into two 
types of the relations; the dispersion relation of the sigma-like mode and that 
of pion-like modes. From the dispersion relation of the sigma-like mode, a 
fluctuation in the direction of the condensate, we obtain only meson-like ex- 
citations with mass gap, but no collective phonon mode. The effect of the 
anomalous Wigner functions, g ~ (aa) , plays an important role for preventing 
the solution of the dispersion relation from the violation of causality as in 
the case of one-component theory. On the other hand, the dispersion relation 
of the pion-like modes, for fluctuations in the direction perpendicular to the 
condensate, contains massless collective mode which can be interpreted as a 
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missing Nambu-Goldstone mode. 

We examined how these mesonic collective modes changes their characters 
in the absence of the exact chiral symmetry. It was shown that the Nambu- 
Goldstone mode becomes massive and is non- analytically transformed through 
kink to the pion mode in the symmetric phase. Such a kink always appears at 
the threshold of continuum of the underlying quasi-particle excitations. 

Finally, we like to make a few comments on the improvement of our for- 
malism and its applications to more realistic physical situations. We have 
derived quantum kinetic equations without collision terms using Gaussian 
Ansatz for the initial density matrix. This approximation corresponds to a 
neglect of all correlations in the system. Justification of this approximation 
remains an important problem for further study. Recently, a modified self- 
consistent Hartree approximation in finite temperature scalar field theory 
has been proposed in the context of chiral phase trantsion[28] and Bose- 
Einstin condensation [29]. Although quasi-particle continuum may be modified 
by these more sophisticated version of the Hartree approximation we expect 
some features may remain unchanged by such improvements. The collision 
terms is needed to describe relaxation phenomena accompanied by entropy 
production[30, 31, 32, 33]. Random rescattering of pions may still play some 
role in the space-time evolution of the chiral condensate and the freeze-out 
process. In this regard, it is interesting to see how our "coUisionlcss" pion 
mode is converted to the Son-Stephanov pion mode [34] which is an isospin 
counterpart of the hydrodynamic spin wave[35]. Our present formalism em- 
phasizes the role of the mesonic mean field in the final state interaction. We 
like to apply our formalism to the kinetic freeze-out dynamics in consideration 
for the effect of the expansion of the system and solve the space-time evolution 
of the condensate assuming boost invariance of the system. It is interesting 
to see whether the Vlasov term might generate strong flow effect due to the 
gradient of pionic mean fleld. This question may be studied by the method 
developed in this work. 

One of the important features of our formalism is that it can describe the 
chiral phase transition both in equilibrium and far out of equilibrium, including 
freeze-out process into free streaming particles, in a unified fashion. This is 

desirable for extracting observable consequence of the chiral phase transition 
to confront with experiments of high energy heavy ion collisions which are 
very complex non-equilibrium phenomenon. 
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A General relations of the Wigner functions 



General relations between the four components of the Wigner functions which 
follow immediately from the definitions (2.21), (2.22), (2.23) and (2.24) are: 



i=;*(p,k,t) = F,,(p,-k,t), (A.i) 

i^;.(p,k,t) = F,,(p,-k,t), (A.2) 
(5,,(p,k,t) = G^(p,-k,t), (A.3) 

where the asterisk (*) stands for the complex conjugate. From the equal-time 
commutation relations of the quantum fields, (2.19) and (2.20), we obtain the 
additional relations between the Wigner functions 



Fij{p, k, t) = Fji{-p, k, t) + SijSk,o (A.4) 

G,,(p,k,t) = G,,(-p.kt) (A.5) 

Gij{p,k,t)^Gji{-p,k,t) (A.6) 

We define the Fourier transforms of the Wigner functions as 



(A.7) 
(A.8) 
(A.9) 
(A.IO) 



Then (A.I) and (A.2) imply that the diagonal components /jj(p, r, i) and 
/ii(p, r, i) are both real functions and (A.4) imply that 



fij{P,r,t) = fji{-p,r,t) + Sij (A.ll) 

while (A.3) implies that gij{p,r,t) and gji{p,r,t) are complex conjugate to 
each other: 



g*^{p,r,t)^gji{p,r,t) (A.12) 



fij{P,r,t) 


k 


-^^•^F,,(p,k,t), 


fij{P,r,t) 


k 


-'"^■^FijipXt), 


9ij{p,r,t) 


k 


-•'•^G,,(p,k,i), 




=Ee- 

k 


-'^•■-G,,(p,k,t), 
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B Equation of motion of the Wigner functions 



The equation of motion of each component of the Wigner functions is given 
below: 



d 

^^^ii(p, k, t) = -K,p+| - a;^. p_k)Fij(p, k, t) 

F,,(p + |, k + q,t) + Q,(p + |, k + q,t) 



, F,;(p-f , k + q,t)+G',;(p-f , k + q,t) ^^ ^ 

^ . ^^%,q 



(B.l) 

d 

Q^Gijip, k, t) = (a;, + o;^. p_k)G,j(p, k, t) 

Gij{p + I, k + q, t) + Fij(p + I, k + q, t)) 



+ E An 



jZ,q" 



G.Kp-I, k + q,t) + F,;(p-3, k + q,t) 



An 



d - 

?^Gij(p, k, t) = -Kp+i + p_k)G'»j(p, k, t) 



+ E -Ana,q 



(B.2) 

(^^.•(p + f , k + q, t) + Fijip + |, k + q, t) 



G'iKp-q/2, k + q,t) + F,/(p-q/2, k + q,t) 



An 



^ 2uji, p-k/2 ^ 2a;j, p-k/2-q 

(B.3) 



a - 

^^^ii(P' k, = Kp+I - p_k)Fij (p, k, t) 



^Y^LyI, (^'^(P + 1, k + q, t) + g,,(p + f , k + q, t)) 
j'A '^'^ V'2a;i,p+t^2u;z,p+|+q 

F,;(p - q/2 , k + q . + G,/(p - q/2 , k + q , t) ^ ^ 1 
^2u;i,p_k/2^2u;j,p_k/2-q J 

(B.4) 
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C Kinetic equations of 0{N) model in the long wavelength approx- 
imation 



Here we present the explicit form of the equations of motion of the Wigner 
functions in the long wavelength approximation: 



1(1 



2 ycoi^p ojj^p^ 



E 



1 / 1 



4 V^'p 4,P. 



1 f 1 



4 V^'p 4,p. 



Uif{r,t)p ■ (iVr/i'j(p,r,t)) 



+ -{-tV,Uiy{r, t) ■ Vp/,v,(p, r, t)) 



-\Uif{r,t)gj'j{p,r,t) - ] [ -i 



4 \^i,p ^f,p. 



P • (iVr%(r,0)^j'j(p,r,t) 



1 / 1 



+ 



4 V^'p 4,p. 



Uij'{r,t)p- {i^r9j'j{p, r, t)) 



+\ Uj'j{r, t)fif{p, r,t) + ^l—^ ^ ) p • (iVrC//j(r, t))/i/(p, r, t) 

+7 ( 4- + -4-] Ufj{r,t)p-{iVrfij'{p,r,t)) 

+ 1 t)gij'{p, r,t) + ] - ) p • (iVrC//j(r, i))^y'(p, r, 



4 V^/,p ^|',p. 



+ 



1 / 1 



+ 



4 V^'p 4,p. 



Ufj{r,t)p- (iVrfifij/(p,r,t)) 



--(-iVrt//j(r, • VpQifip, r, i)) 



(C.l) 
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ifijiP, r, t) = +{uji,p - ujj,p)fij{p, r, t) + 

-E 



-\Uij'{r,t)fj>j{p,r,t) -] 



4 V^'p 4,P. 



I ( — + ^)p-(^Vr^(p,r,t)) 



p- {iVrUij'{r,t))ffj{p, r,t) 



1 / 1 



+ 



4 V^'p 4,p. 



Uij'{r,t)p ■ (iVr/i'j(p,r,t)) 



+ -(-iVr^i,'(r, t) • Vp/,v,(p, r, t)) 



-\Uij'{r,t)gfj{p,r,t) -] 



4 V^'p 4,P. 



P ■ {i^rUij'{r,t))gfj{p,r,t) 



1 / 1 



+ 



4 V^'p 4,p. 



C/ijv(r,i)p • (iVr5j/j(p,r, i)) 



+ 2 (-^VrC/i/(r, • VpQfjip, r, t)) 



f - 1/1 

+1 t^j'j(r,*)/i/(P,r,t) + - 



4 V^'p '^i',?. 



p- (iVrC/,'j(r,t))/ij/(p,r,i) 



+ 



1/1 



+ 



4 V^'p ^hp. 



Uj'j{r,t)p- (iVr/ij'(p,r,t)) 



--(-iVrC//j(r, • Vpfifip, r, i)) 



+ 1 Uj'j{r,t)gij,{p,r,t) + 



III 



+ - 



1/1 



+ 



4 Vo;, 
1 



p- {iVrUj'j{r,t))gijr{p,r,t) 



4 V'^'p ^},p. 



Uj'j{r,t)p- {iVrgij'{p,r,t)) 



-^{-i^rUfj{r, t) ■ VpC/yv(p, r, t)) 



(C.2) 



44 



Wij{p, r, t) = -{uJi,p + ujj,p)fij{p, r, t) - 

+E 



1/1 



2 \a;i,p a;j,p^ 



P ■ (^Vr5ij(p,r,t)) 



1 / 1 

%(r,t)5i'i(p,r,t) - - 



4 V^'p 4,p. 



P ■ (^Vr^7ii'(r,t))^j/j(p,r,t) 



1 / 1 



+ 



4 V^'p 4,p. 



^ii'(r,t)p • (iVrfifj'j(p,r,t)) 



+ -{-iVrUif{r,t) ■ Vpgj'j{p,r,t)) 



+ \Uij'{r,t)fj,j{p,r,t) -] 



1 / 1 



+ 



4 V^r,p ^j',p 
1 



P ■ (iVr^7j/(r,t))/,/j(p,r,t) 



4 V^'p 4,p. 
1 



C/ijv(r,i)p • (iVr/j'j(p,r,i)) 



+ 2(-^VrC/i/(r,t)-Vp/,.,-(p,r,i)) 



f 1/1 
+ 1 Uj>j{r,t)gij>{p,r,t) + - 



1 



4 V^'p 4,P. 



P- («VrC/j'j(r,i))^yv(p,r,t) 



+ 



1/1 



+ 



4 V^'p ^hp. 



Uj'j{r,t)p- {'i^r9ij'ip,r,t)) 



--(-iVrC//j(r, • Vpc/yv(p, r, i)) 



+ 1 Uj'j{r,t)fijr{p,r,t) + 



III 



4 V'^'p 4,p. 



p- (iVrt/,'j(r,i))/ij/(p,r,i) 



+- 



1/1 



+ 



4 V'^'p ^',p. 



[/j/j(r,i)p- {iVrfij'{p,r,t)) 



-^{-i^rUfj{r, t) ■ Vpfij'ip, r, t)) 



(C.3) 
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Wij{p, r, t) = +(wj,p + ujj,p)fij{p, r, t) + 

-E 



1 / 1 

%(r,t)c/i'i(p,r,t) - - 



4 Vo;?. 





1 






1 






1 0- 



P- («Vr5ij(p,r,t)) 



p- («Vr^/ij'(r,i))^j/j(p,r,t) 



1 / 1 



+ 



4 V^'p 4,p. 



+ -(-?Vr^7i/(r,t) • Vpgj'j{p,r,t)) 



+ \Uij'{r,t)fj,j{p,r,t) -] 



1 / 1 



+ 



4 V^r,p ^i',p 
1 



P ■ {i^rUij'{r,t))fj'j{p,r,t) 



4 V^'p 4,p. 
1 



C/ijv(r,i)p • (iVr/j'j(p,r,i)) 



r _ 1/1 

+ 1 Uj>j{r,t)gij>{p,r,t) + - 



1 



4 V^'p 4,P. 



P- («VrC/j'j(r,i))^yv(p,r,t) 



+ 



1/1 



+ 



4 V^'p ^/,p. 



--(-iVrC//j(r, • Vpc/yv(p, r, i)) 



+1 t^j'j(r,i)/y'(P>i',i) + 



1/1 



4 V'^'p 4,p. 



p- {iWrUjrj{r,t))fij>{p,r,t) 



+ - 



1/1 



+ 



4 V'^'p ^',p. 



[/j/j(r,i)p- (iVr/ij'(P>r,i)) 



-^{-i^rUfj{r, t) ■ Vpfij'ip, r, t)) 



(C.4) 



D Linearized Vlasov equations of 0{N) model 



The complete set of the linearized Vlasov equations discussed in section 3 is 
listed below: 
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i6fij{p, r, t) = -{ui^p - LJj,p)6f^j{p, r, t) - + — )p • ^rSfij 



+ (/eq.Kp) -/eq.K.p)) 



5C/,,(r,i)-l(-^--L)p 



+ 2 Vp(/eq.(^i,p) + /eq.K,p)) ' ^rSUij 



— (U!. 



-'^i,p)^/ii(P,r,t) - ^(-^ + — )p 



2 ^a;j_p a;j_p 



-(/eq.(^i,p) - /eq.(^i,p)) 



4^<P ^Ip 



:Vp(/eq.(a;i,p) + /eq.(a;j,p)) • VrSU^ 



iSgijip, r, t) = {uji^p + a;j-p)%j(p, r, t) + ^-)p • VrSgij 

z ^oji^p ujj^p' 



+ (/eq.Kp) +/eq.K-,p)) 



4 ^Ip 



+ 2Vp(/eq.(^j,p) - /eq.(^j,p)) ' Vr^C/y 



i / 1 

«%(P, r, t) = -{uji^p + ujj^p)6gij{p, r,t) - - I 



2 ^uji^p ujj^p 



P ■ ^r^g 



+ (/eq.Kp) +/eq.(^i,p)) 



6U,,{r,t) + U-^--^)p 



4^<P ^Ip 



-^'i^p{feq.{^i,p) - /eq.(^i,p)) " ^rSUij 



The Fourier transform of these equations yields: 
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{ 



UJ -\- UJi — UJj — 



p • k/ 1 1 



(- + -)W«,(p,k, 



00 



{l + ^(^ - ^)P ■ k}(/eq.(0 - /eq.(j)) " • Vp(/eq.(0 + /eq.(i)) 



(D.5) 



\uj - UJi + UJj + ^77^ (— + —)\sfij{p,k,u) 

I Z ^UJi UJj' ) 



5t/ij(k,a;) 
(D.6) 



(uj- UJi -UJj + ^^^(— - —)\6gij{p,k,uj) 



p ■ 1 1 

p • k|(/eq.(0 + /eq.(j)) ' ^k ■ Vp(/eq.(0 " /eq.(j)) 







^ 1 













5t/ij(k,u;) 
(D.7) 



+ + + ^-^(^ -)\5gij{p,'k,uj) 

I Z ^UJi UJj ' ) 



(5C/y(k, a;) 
(D.8) 



where we have used shghtly abbreviate notations: uJi — uji^p, ujj — ujj^p, /eq.(0 

/eq.('^i,p) and /eq.(i) = feq.{^j,p)- 
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